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We propose a scheme to probe the non-Abelian statistics of the collective anyonic excitation in Kitaev's 
honeycomb model with cold atoms in an optical lattice. The generation of the anyonic excitation can be realized 
by simple rotating operation acting on an effective spin-1/2 system, which is encoded in the atomic hypcrfinc 
energy levels. The non-Abelian nature of the anyonic excitation is manifested by the braiding of four vortices, 
which leads to different operations on the subspace of degenerate ground states, and thus results in different 
final states. Here, by introducing an ancilla atom, the effective control over the lattice atoms can be realized 
and the final different states can also be imprinted on the ancilla and further distinguished by measurement. 
© 2013 Optical Society of America 
OCIS codes: 000.1600, 140.3325, 270.5585 



Quantum statistics is one of the basic concepts in quan- 
tum physics [1]. It's also the root of some condensed 
matter physical phenomena, such as, superfluid, Bosc- 
Einstein condensates, etc. In three-dimension, due to the 
rotating and time reversal symmetries, particles obey ei- 
ther the Bose-Einstein or Fermi-Dirac statistics corre- 
sponding to their statistical phase a (a = 0/ir) to be or 
1. But, there exist cases with a ^ or 1 corresponding to 
more exotic particles called anyons, which maybe cither 
Abclian or non-Abelian because of symmetry breaking 
in 2-dimcnsion [2] . The elementary excitation with frac- 
tional charge in fractional quantum Hall effects (FQHE) 
is anyon satisfying exotic statistics [3,4], which provides 
a good platform to understand anyonic statistics. For 
Abelian anyons, the anyonic wave-function acquires a 
phase factor coming from the winding of two quasiparti- 
cles, which is just the one-dimensional representation of 
the braid group Bn (for A*" indistinguishable particles). 
For non-Abelian anyons, quasiparticlcs exchange would 
be represented by a matrix acting on the subspace of 
the degenerate ground states, which is an exactly non- 
commutative representation of the braid group in higher 
dimension [5]. In particular, Kitaev showed [6] that non- 
Abelian anyon is a promising candidate for topologi- 
cal quantum computation, which is fault-tolerant at the 
hardware level. 

Up to now, no experiment has unambiguously verified 
the existence of Abclian or non-Abelian anyons although 
promising progress have been made in v = 5/2 FQHE 
state [7] . This is mainly due to the extremely complicated 
noise channels in a quantum many-body solid-state sys- 
tem. On the other hand, a cold atom system possesses 
unprecedented possibility of controlling almost all rele- 
vant physical parameters [8] and thus been recognized as 
an ideal system for quantum simulation. Proposals of ob- 
serving the non-Abelian statistics are also presented with 
cold atoms in s-wave [9,10] and p-wave [11] superfluidity. 



For Kitaev's toric model, it has been proposed to simu- 
late the anyonic interferometry in virous systems [12-19] 
with experimental verification in optical systems [20-22] 
using a method of generating dynamically the ground 
states and the excitations of the model Hamiltonian. 
However, there is few reports about the verification of 
non-Abelian statistics in the Kitaev honeycomb lattice 
model [23]. Recently, theoretical implementation of this 
model has been proposed with cold atoms in an honey- 
comb lattice [24]. Meanwhile, elementary experimental 
verification of the Dirac physics associated has also been 
made with cold atoms in a honeycomb lattice [25] based 
on theoretical proposals in [26,27]. 

Here, motivated by the above advances, we propose 
a concrete example of probing the exotic non-Abelian 
statistics in the Kitaev honeycomb model with cold 
atoms in optical lattice. The honeycomb lattice is con- 
structed from cold atoms located in an optical lattice 
and the anisotropic spin interaction is induced from spin- 
dependent tunneling between the nearest neighbor of 
the lattice [24]. To get effective control over the lat- 
tice atoms, following Ref. [14], we introduce an auxiliary 
atom. Then, creation and manipulation of anyons are 
implemented by Rydberg gates based on dipole-dipole 
interaction between the ancilla and lattice atoms. Fi- 
nally, the exotic non-Abelian statistics can be recorded 
by the ancilla during the implemented braiding process. 
Therefore, detect the final state of the ancilla can reveal 
the non-Abelian statistics. Furthermore, it is usually dif- 
ficult to read out the states of the anyons due to the fact 
that they are degenerated and neutral. Therefore, we in- 
troduce an auxiliary atom to imprint the difference of the 
anyonic final states after braiding. Meanwhile, for prob- 
ing the non-Abelian statistics, one needs at least four 
anyons. If only two anyons are involved, the final state 
after braiding will be different from the initial state by 
a phase factor. However, this difference is insufficient to 
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Fig. 1. Implementation of the Kitaev model, (a) The 
model is a honeycomb lattice model of spin 1/2 system 
in a plane. The index of the 6 atom within a plaquettc is 
indicated, (b) There are 3 different interaction between 
the nearest neighbor of the lattice along the x — x, y — 
y and z — z directions, (c) The pseudo spin states are 
encoded by | \) and | \) of an atom, an excited state |e) 
is also introduced. Two coupling laser beams L\ and L 2 
are used to realize the interaction in three independent 
directions as we needed. 



verify that the involved anyons are of non- Abclian nature 
since braiding of Abelian anyons may also cause such a 
phase difference. 

Let us begin with a brief review of the Kitaev hon- 
eycomb lattice model, which is an anisotropic spin-1/2 
model with different nearest neighbor interactions in dif- 
ferent bonds [23] 



H = - J * E °*°i- J v E °K- J * 

y— links 



x — links 



E rfrf 

z — links 



>i"k, 



. (1) 

where a" with v G {x, y, z} being the pauli matrix on 
the site j, J„ are coupling strength for different links as 
indicated in Fig. 1. To get the ground states, diagonal- 
ization of the model H is needed, which can be greatly 
simplified by constructing a set of integrals of motion 
W p = o\<j\o\o%<j\o% with p being a label of a plaque- 
tte of the lattice and the spin index is labeled as in Fig. 
1(a). Expanding every spin into a four Majorana oper- 
ators notation, the model Hamiltonian in the enlarged 
space reads [23] 
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where Ajk = 1J a Ujk if spins j and k are connected, 
otherwise Ajk = 0. Here, Ujk = ib^b^ with bj and Cj 



being Majorana operators, [H, Ujk\ — 0, uL 



Ujk, 
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h ik 



1. We first divide the total original Hilbert space 
into different sector which is eigenspace of W p . So W p 
can be transformed to the eigenspaces of a set of states 
given by a certain configuration of eigenvalue of Ujk- The 
connection between w p (eigenvalue of the W p ) and Ujk 
(eigenvalue of u jk ) is w p = U(j,k)&oundary(p) u ok- The 
variables Ujk and the numbers w p can be interpreted as a 



Z2 gauge field and the magnetic flux through the plaquc- 
tte p, respectively. We say that the plaquette p carries a 
vortex when w p = — 1 and the ground states are denoted 
by the vortex- free sector, i.e. all w p = 1. Then go to the 
momentum representation via the Fourier transforma- 
tion, we can get the ground energy, ground states and the 
phase diagram. In the triangle region of | J x | < | J y | + 1 J z \, 
\J y \ < \J Z \ + \J X \ and \J Z \ < \J X \ + \J V \, it's gapless phase 
B with non- Abclian anyonic excitation [23]; otherwise, 
it's gapped phase A with Abelian anyonic excitation. A 

magnetic field V = — Ylj(hx a j + ^2/ cr J + ^z°j) ^ s intro- 
duced to open a gap between the ground and the ex- 
cited states in phase B [23]. This gap will exponentially 
suppress the interaction between vortices and fermionic 
modes near the singularity of the spectrum, and thus 
makes the statistics of non- Abelian anyon well-defined. 
Now, we turn to construct the honeycomb model of Ki- 
taev using ultracold fermionic atoms in an optical lattice. 
Following Ref. [24] , we consider a cloud of ultracold Rb 8 ' 
atoms trapped in an optical lattice which are formed 
by interfering standing laser beams. First, completely 
suppress the tunneling and the spin interchange in the 
vertical Z direction by raising the corresponding barrier 
potential, and thus the atoms distribution form a two- 
dimensional configuration in the X — Y plane. Assuming 
that the system in the Mott insulator regime and the 
atomic density is roughly one atom per site. Each atom 
is treated as a two-level system, labeled with the effective 
spin index a =t, I, with an energy splitting of 6834 MHz. 
To form the honeycomb lattice, we apply trapping poten- 
tial Vj (x, y) = Vosin 2 [k» {x cos 6j + y sin 9j ) + 4>q] formed 
by a pair of blue-detuned traveling laser beams above the 
X-Y plane with an angle </>ii, where, j = 1, 2, 3, #1(2,3) = 
7r/6(7r/2, -7I-/6), 0|| = 2arcsin(l/V3), and <rj = tt/2. To 
engineer the wanted interaction in different directions, 
applying the following three blue-detuned standing- wave 
lase beams V ua {x, y) = V va sin [fc(xcos6^ +ysin6^)] in 
the X-Y plane along certain directions, where v denote 
the tunneling direction x-x, y-y, z-z, 9' x _ x ~ — w, 0' y _ y = 
ir, 9' z _ z = 7r/3. These laser beams will introduce spin- 
dependent potentials V va = V„ + |+), y (+| +K,_|— }„(— |, 
where |+)„ and |— } v are the eigenvectors of the eigenval- 
ues 1 or —1 of Pauli operator a v (cf. Fig. 1), respectively 
That is to say, different spin directions will feel differ- 
ent effective potential, and thus have different tunneling 
rates. The potential depth between the nearest neigh- 
bor along the v direction can be changed from Vb/4 to 
Vb/4 + V va [24]. The parameters V v ± can be tuned by 
varying the laser intensity of L\ and L 2 in the v direc- 
tion. Here, we are interested in the regime t u+ 3> £„_, 
i.e., the atom can virtually tunnel only when it is in the 
eigenstate of |+)„. Then, the spin-dependent potentials 
will lead to anisotropic Ising exchange interaction for 
each tunneling direction, and the total interaction is in 
the form of [24] 
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Fig. 2. (color online) The creation and manipulation of non-Abelian anyons. (a) Initial state of the ancilla atom (the 
red dot) is |+) a ex |0) Q + |l) a and the solid blue circle denote a vortex excitation bound to a plaquette. When moving 
the ancilla close enough to the code atom, 2-qubit U z gate is achieved, (b) The manipulation of the anyons, only 
when the ancilla is in state |l) a , an effective operation will act on the atom transforming the w p from one to the 
other sector, and thus the vortex configuration can be changed as we needed. 



In the case of K-/K+ > 1 and U+ « f_ ra £/+_ w U, 

where U a and U-\ are the intra- and interspin coupling 

energy, J v = t 2 +u /{2U) and h v = 4t„ + /U. Then, we get 
a Hamiltonian of Eq. (1) with additional effective mag- 
netic field terms hv^.a^, which can be used to open 
the gap in phase B. Therefore, this implementation nat- 
urally leads to the effective magnetic field terms, and 
thus makes it more easier for realization. 

To probe the non-Abelian statistics, we need at least 
two pairs of non-Abelian anyons due to their supcrsc- 
lcction rule. The generation of a pair of anyons can be 
realized by a simple a z rotation acting on a single spin, 
which is in a vortex- free configuration. But, the gener- 
ation process will accompany with fermionic excitations 
in the mean time. Fortunately, the unwanted fermionic 
excitations can be exponentially suppressed by the en- 
ergy gap opened by the effective magnetic field. The a z 
rotation may be implemented by laser-atom interaction, 
but single atom addressing is difficult experimentally. To 
avoid such difficulty [14], we consider an ancilla atom 
that is trapped above X — Y plane by three standing 
wave laser beams along three direction with its inter- 
nal Zeeman levels |0) a and |l) a . The a z rotation can be 
realized by Rydberg gates based on dipole-dipole inter- 
actions between the ancilla and a lattice atom [28,29], 
which maximally reduce the addressing difficulty and 
eliminates the need of cooling both the lattice and an- 
cilla to the physical ground states in the cold collisions 
implementation [30]. On the other hand, this single an- 
cilla atom can be moved by trap potential without de- 
coherence [31] and this is important to guarantee our 
implementation. As a result, 2-qubit unitary operations 
Uf = |0) a (0|®/+|l) a (l|®of acting on site atom i and 
the ancilla atom are implemented, which are sufficient 
to create and manipulate the anyons [23], as shown in 
Fig. 2. Experimentally, an arbitrary qubit hopping time 
is 8/is with a nearly unit fidelity [32]. On the condition 
that the ancilla atom is in the |l) a state, erf rotations 



can be implemented on site i atom. The generation of 
two vortex configuration can be represented as 



H 2 - v = vfHaf 
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where a z denotes a <j z rotation on site i, (i,j) label a x- 
link, and (i, k) labels a y-link. meanwhile, we can also get 
two vortices via a of or af rotation but with different 
anyons configuration. This simple configuration due to 
Uu = — 1 at the site i along the z — z direction with the 
neighbor honeycomb lattice u p — — 1, and the vortex- 
free configuration as it gives lo p — 1 for all plaquettes p. 
The exotic collective excitation are non-Abelian anyons 
bound to the vortices, and the non-Abelian anyons are 
Majorana fcrmions [23]. At low temperature, the initial 
state of the pair of Majorana fermions is typically a vac- 
uum state |0) [11]. 

To make our scheme more robust to decoherence and 
avoid the problem of single atom addressing, we use 
the conditional generation of anyons using U". When 
a = z, the ancilla atom, vortex- free ground state \gs) 
and 2-vortex state |0) are entangled as shown in Fig. 2, 
i.e. |0) a |c/s) + |l) a |0). Therefore, the exotic non-Ablcian 
statistics can be recorded by the ancilla atom, i.e., we 
can also use the state of the ancilla atom to manifest 
the non-Abelian statistics of anyon. The state of ancilla 
atom can be addressed independently, and thus can be 
measured by standard probing technique. 

We next consider the simplest vortices configuration 
that can experimentally reveal the non-Abelian statistics 
as shown in Fig. 3, which consists of two pairs of vortices. 
To reveal the non-Abelian statistics, we would like to 
move one anyon (Majorana fermion 3) around the other 
(Majorana fermion 2), where the two come from differ- 
ent Majorana fermion pairs. Indeed, this braiding can be 
implemented by the operator s 2 3 = (J t (J l a ?, a 2 a i (J 6 P3]- 
For a system with 4 vortices, the degeneracy of the 
ground states is 4. For an adiabatic braiding of vortices 
2 and 3 once, the effect is a unitary operation acting 
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Fig. 3. (color online) The braiding of two non-Abelian anyons. In a topological sense, a counterclockwise braiding 
(or exchange) of two anyons is equivalent to move a anyon counterclockwise by 2-7T around the other. We choose the 
braiding take place between anyons 2 and 3 due to superselection rule. The operator S23 = erf cfcf erf erf af correspond 
to a looping trajectory which is exactly the loop when anyon 3 move around anyon 2. 



on the the subspace of the ground states, which can be 
represented as [33] 



R * = T* 



1 
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where R{ correspond to anticlockwise exchange between 
two Majorana fermions % and i + 1, and R^ 1 represent a 
clockwise exchange. 

For our particular implementation in Kitaev's honey- 
comb model, we prepare the initial state of the ancilla 
atom in |+) a oc |0) a + |l) a . We first create the lattice 
atoms in the 2-vortex state ji^), i.e., with anyons 1 and 
2. Secondly, bringing the ancilla close enough to the site 
"6" atom so that their collision interaction is in effect, ie., 
implement a U z gate to the atom and the ancilla. There- 
fore, the entangled state |0) a |t>2) + |l)a|^4) is produced 
with 1 1*4) being the 4- vortex configuration state. By do- 
ing so, we have created anyons 3 and 4 on the condition 
of the ancilla in the state |l) a . Because the braiding will 
be implemented between anyons 2 and 3 and the anyons 
1 and 4 are not surrounded in the braiding pathway, so 
it is convenient to work in the subspace of vortex 2 and 
3 on which operation R2 act. Thirdly, move the ancilla 
close to the target atom " 6" again to realize Ug gate op- 
eration, and then "1" to "5" successively to realize the 
braiding operation S23 = C/5 C/|C/fC/|t/^t/g when the an- 
cilla atom is in the |1) state (cf. Fig. 3 ). After braiding 
of anyons 2 and 3 twice, we get |0) a |^2) + e~ 1 ^ \^) aR^i) 
with i^lwi) is a degenerate ground state with a fcrmionic 
occupation. The vortex pair state of 2 and 3 changes from 
the vacuum to a fermionic state with a fermion occupa- 
tion. Finally, repeat the second and first steps to bring 
the system back to the vortex-free state, the final state 
of the system is (|0) a + e~ l ^|l) a )|c/s) with the ancilla in 
the state of |/) oc |0) a + e _l ^|l) a . Then, compare the 



initial and the final states of the ancilla one can find 
that a phase factor — ir/2 is generated due to the non- 
Abelian nature of anyons. The recorded interferometry 
result of the non-Abelian statistics can be read from the 
ancilla atom by a local projective measurement after a 
7r/4 phase gate. The single atom gate and measurement 
can be realized with very high efficiency. 

Note that while the anyon obeys Abelian statistics, 
our procedure will bring the final state of the ancilla to 
|0) — |1) (braiding once); for bosonic or fermionic statis- 
tics, the ancilla's state will stay unchanged. Therefore, 
probe the final state of the ancilla will verify the no- 
Abelian statistics. On the other hand, for no-Abelian 
anyons, when the braiding is repeated twice and four 
times, the final states of the ancilla atom will be the 
same as the Abelian anyon case (braiding once) and be 
identical to the initial state (up to an overall n phase as 
i?2 = —1), respectively. For Abelian anyon, two succes- 
sive identical braiding will lead the ancilla's final state to 
be the same as its initial state. This can serve as further 
distinction of the Abelian and no- Abelian anyons in this 
particular model. 

In summary, we propose a scheme to probe the non- 
Abelian statistics of Majorana fermions bounding to the 
vortex excitation in the Kitaev honeycomb lattice model 
with effective magnetic field using cold atom in an optical 
lattice. The exotic non-Abelian statistics can be recorded 
and detected by the internal state of the ancilla atom, 
the measurement of which can be realized with very high 
efficiency by standard techniques. 
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